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We study the absorption and scattering of massless scalar waves propagating in spherically sym- 
metric spacetimes with dynamical cosmological constant both in low-energy and high-energy zones. 
In the former low-energy regime, we solve analytically the Regge-Wheeler wave equation and ob- 
tain an analytic absorption probability expression which varies with My/A, where M is the central 
mass and A is cosmological constant. The low-energy absorption probability, which is in the range 
of [0,0.986701], increases monotonically with increase in A. In the latter high-energy regime, the 
scalar particles adopt their geometric optics limit value. The trajectory equation with effective 
potential emerges and the analytic high-energy greybody factor, which is relevant with the area of 
classically accessible regime, also increases monotonically with increase in A, as long A is less than 
or of the order of 10 4 . In this high-energy case, the null cosmological constant result reduces to the 
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The absorption and scattering of particles outside a black hole has been a topic of research since 1970s. This 
issue refers to many key theoretical physics problems, such as quantum entanglement Q-!^, entropy 0,0], information 

ravity where fundamental issues about gravity physics and the 



[tHy|J which arise in black hole physics and quantu m g 
structure of space are of interest. After Hawking [lO - 



structure of space are of interest. After Hawking [10h12| ] proved black hole can emit thermal radiation via quantum 
field theory, the absorption and scattering issue is focused gradually in black hole physics. Many people, such as 
Unruh 13], Sanchez 14-16], Das [ljj, Emparanflcl [ijj] etc., have done a lot of works on the black hole scattering in 
the aspects of calculation of absorption greybody factor and cross section as well as black hole's based properties such 
as unitarity and reciprocity T4], phase shifts and diffraction pattern flfjl ]. angular distribution 16], extra dimensional 



effect 



17 



19] and so on. What this paper 



131 ] refers to is the Unruh's work in which the original analytical method to 
calculating black hole cross section is shown in detail in 1976. The radiation equation, i.e. Regge-Wheeler equation 
in perturbation theory, is broken down into three sub-equations in following three regions: one is near the horizon, 
one is intermediate region and the last is far from the black hole. Then these sub-equations can be solved easily by 
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means of well-known functions or specific functions. Finally, the absorption coefficients are obtained by matching 
these solutions in the areas of overlap between the regions. It is found that the absorption cross section of Dirac 
particles is exactly 1/8 of the massive scalar particles. The total cross section of massive scalar field is a function of 
ingoing particles frequency. Scalar cross section a s and transmission coefficients (T s )i 13[ are 
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where v = yx — m 2 jui 2 is the velocity of particle, I is the angular quantum number, m is the mass, to is the energy. 
Since then, this analytical method has become a standard procedure to calculate absorption cross section. 

Meanwhile, there is also the others method to obtain the cross section, for example in 1978 an algebraic method, 
i.e. power-series approximate expansion, is improved by Sanchezv 



14 



Anderson developed a phase-integral method 



16] in Schwarzschild space. Additionally, 
20j in 1995. On the aspect of modeling methodology, a wide variety of 



scattering models have been proposed including fermion emission 
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Neutrinos 



charged leptons 



24J, elastic scattering 



25j, D-brane spectroscopy 26j [2jj and so on. Attentively, in 1997, Das, Gibbons and Mathur [17| 



calculated the low energy absorption cross section of massless scalars and spin 1/2 particles which immerses a general 
asymptotically flat black hole with arbitrary dimensions. The former cross section is equal to the area of the black 
hole and the latter one is related to the area of true metric. Unlike this work, a non-asymptotically flat de Sitter 
universe is considered here. 

Recently, more and more attentions paid to high dimensional black hole such as high dimensional Schwarzschild-like 
black hole with analytical 
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and numerical 
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3l| methods, high dimensional rotating Kerr-like black holes 



high dimensional Schwarzschild-de Sitter-like black hole [35j and so on. Specially, in the last work [35j the 



absorption and scattering of multi-horizons space was first derived with analytical and numerical techniques by Kanti, 
Grain and Barrau in the view of high dimensions . It is shown that Hawking radiation can reveal much valuable 
information such as the space dimensionality, the space curvature, the bulk cosmological constant. Here, we also 
consider a multi- horizons space but in the view of usual 4 dimensions by the similar method in Ref. js5|. 

On the other hand, the accelerating expanded universe is justified by more and more astronomical observations 



such as la Supernovae (SNe la) 38J-42], cosmic microwave background (CMB) 
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43 461 in the Wilkinson Microwave 



Anisotropy Probe (WMAP). The expanded universe shows us a nonzero cosmological constant A 
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and the 



topology structure presents a de Sitter one. Like the standard null cosmological constant gravity, a lot of literatures 
focuses on Hawking radiation in an inflationary universe including two-dimensional toy models 49h51| , radiative tails 



[52j, degenerate special Nariai black hole j53] and so on. In these works, Brady, Chambers, Krivan and Laguna [52 1 
demonstrated the existence of exponentially decaying tails at late times outside SdS and Reissner Nordstrom-de Sitter 
(RNdS) black holes where the S wave mode asymptotes to a nonzero value. In this work, an advanced time v = t + r* 
and a retarded time u — t — r* where r* is the tortoise coordinate is used. However, there is no work given out 
the explicit expressions of analytic absorption probability with cosmological constant as far as we know. The reason 



The usual 4D multi-horizons space scattering was firstly studied by Giirsel, Sandberg, Novikov, and Starobinsky with charged Reissner 
Nordstrom black hole [36l,l37ll. 
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maybe is the complex ordinary differential equation which could not be solved completely and do not give us the 
general solution. Considering these situations, we adopt the approximate treatment with low-energy and high-energy 
limits assumptions and attempt to calculate the absorption and scattering cross sections or greybody factors outside 
SdS black hole. So the whole problem is divided into two parts: low-energy regime and high-energy regime. In the 
former case, the similar Unruh's method is used. But we will only use two zones, i.e. near event horizon and 
cosmological horizon respectively, and match the solution in the intermediate zones. In the latter high-energy case, 
we use the innermost stable circular orbit to fitting the high energy particles and try to obtain the greybody factor. 
Another purpose of this paper is to inspect whether the cosmological constant can influence the cross section of usual 
Schwarzschild black hole. 

Finally, we want to clarify two points: one is the units of cosmological constant A and the other is the polarization 
of absorption cross section in black hole. For the former case, we use the Planck Units (G=c=fi,=l) to specify A in 
this paper. In this units, A has dimensionless form. Also, in order to simplify the calculation of graybody factor, 
the mass of black hole M is adopted to unity in this manuscript. Hence, if we convert cosmological constant A to 
ordinary units, we can multiply it quantificationally by its conversion factor via A — > A \GM/c 2 ) . So in the ordinary 
units, A has a dimension of m~ 2 . About these specific details, please see the Refs. 5J-[56(. For the latter polarization 
case, incoming scalar particles have no polarization. In fact, the polarization exists in the scattered-out beam for the 
charged particles outside a black hole, despite the incoming beam is not polarized. To take Reissner-Nordstrom black 
hole for an example, since the effect of absorption, the non-zero scattering amplitude provides a polarization effect. 
Hence, the scattered-out beam is polarized, and the polarization vector is perpendicular to the scatterin g pl ane. The 
polarization grade is the function of azimuth angle. About these specific details, please see the Refs. 0, 13]. 

This paper is organized as follows: in section II, we present the Schwarzschild-de Sitter black hole solution and 
obtain the general emission equation. In section III, we calculate absorption probabilities for scalar particles emission 
in low-energy regime. In section IV, we use the geometric optics limit and calculate the greybody factor of scalar 
particles emission in high-energy regime. Section V is a conclusion. We adopt the signature (+, — , — , — ) and put h, 
c, and G equal to unity. 



II. SCHWARZSCHILD-DE SITTER BLACK HOLE AND GENERAL EMISSION EQUATION 



The static spherically symmetric metric of the Schwarzschild-de Sitter space 

1 



59| is 



ds 2 = f{r)dt 2 



-dr l 



which is obtained via solving the Einstein gravity field equations coupling with cosmological constant A, 



(3) 



(4) 



The central black-hole mass M is contained in metric function f(r) = 1 — 2M/r— Ar 2 /3. In this paper, the cosmological 



constant A is considered as a free parameter as shown in Refs. (60l463| . The interesting features for the space metric 
([3]) are its two bounded horizons, one is inner horizon (i.e. black-hole horizon) and the other is outer horizon (i.e. 
cosmological horizon). Under the limit A — > 0, this metric has exactly the same line-element as the Schwarzschild 
space. But in the limit of M — > 0, it reduces to the de Sitter one. By using these horizons, the metric function can 
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FIG. 1: event horizon r e (dashed line), cosmological horizon r c (dotted line) and interval \r c — r e \ (solid line) versus cosmological 
constant A in the range of [0, 0.11] for M = 1. 

be rewritten mathematically as a new form, 

f(r) = ^~(r - r e )(r c - r)(r - r ). (5) 

The singularity of metric ([3]) is determined by equation f(r) — 0, which is also the null hypersurface condition for 
static spherically symmetric space. So, these solutions are inner horizon r e and outer horizon r c , as well as a negative 
solution r Q = — (r e + r c ). The last one has no physical aspects, and is not considered here. The former two positive 
solutions r e and r c are listed as follows, 

2 2 
r c = —= cost/ and r e — —= cos(27r/3 — rj), (6) 
VA ' VA 

where r\ = 1/3 arccos(— 3MVA), (rj £ [7r/6 < 77 < 7r/3]). The real physical solutions are accepted only if A satisfies 
AM 2 < 1/9 64|. If the cosmological constant A reaches its maximum, the Nariai black- hole will appear, which has 



attracted a lot of academic attentions 



65rH- 

Here we plot the curves of event horizon r e and cosmological horizon 



r c in FigJTJ which show clearly that the outer cosmological horizon r c decreases with increase in cosmological constant 
A. Else, r c decreases exponentially quickly when A is under the order of 10 -2 . And the amplitude of variation for r c 
is relative slower when A is larger than 1CP 2 . But the inner event horizon r e gently increases oppositely with increase 
in A. On the other hand, it is known that our world is located between event horizon and cosmological horizon. So we 
also plot the interval \r c — r e \ versus cosmological constant A in the same FigH] It illustrates that with the increasing 
A, the interval range decreases monotonically. Else, the rate of decreasing is very large both in the primary stage and 
the terminal stage of dynamic cosmological constant. 
Employing the factorized ansatz 

Q = -jL=-R u (r,t)Y lm (0, t l>), (7) 

where Yi LJ (9,(j)) is the spherical harmonics function, the minimal coupling massless scalar field <&(t, r, 6, (j)) equation 

1 d ( .... d 
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is broken into spherical harmonic equation (SHE) and the radial radiation equation (RRE). The former SHE is one 
of the main characteristic for sphere-symmetric space, and is not shown here. The latter RRE contains much of the 
information for Hawking radiation, which is shown as, 



1 



d 2 (R 



d 



2 /(r 



d 



r 



-1(1 + 1) 



Re 



= 0, 



(9) 



f(r) dt 2 \ r J ' dr \ J w <9r 

where i? w (r, i) is the time-dependent radial function. Radial radiation Eq.® is the propagated master equation ((9J, 
which determines the evolution of evaporating black-hole for scalar particles. 

Eliminated the time variable by the Fourier component e" lujt via R u (r,t) — > ^ u i(r)e~ lut , Eq.® can be rewritten 

as 



whose dominion potential is given by 



V(r) = f(r) 



ldfjr) | 1(1 + 1) 
r dr r 2 



(10) 



(11) 



Here in order to inspect the influences of angular momentum and cosmological constant on this potential, we plot 
Figs[2] and [3] respectively. From these figures, we can see easily that the potential V(r) depends on parameters A and 
I sensitively. It is illustrated that potential V(r) is higher with increase for the angular momentum quantum number 
I from Figl21 and it also suggests that for the spherically symmetric s wave, i.e. I = 0, its absorption capacity is the 
strongest than the others. So, the dominant s wave is only considered in this paper. But viewing from FigO the 
potential becomes lower with increase in cosmological constant. 
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FIG. 2: potential of Schwarzschild-de Sitter black hole 
with cosmological constant A = 1CP 4 for unit M — 1 and 
various angular momentum wave. 



FIG. 3: potential of Schwarzschild-de Sitter black hole 
with various cosmological constant for unit M = 1 and 
angular momentum I = 1 (p wave). 



Now we introduce the tortoise coordinate 



r* 



1 

2M 



dr 

Wr 



(12) 
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The tortoise coordinate can be expressed by surface gravity as follows, 



r* = 



1 

2M 



i ln ( 1- ^ 



where the definition of surface gravity is Ki = i 
we can obtain following explicit expressions of corresponding three surface gravities as 

(r c ~ r e )(r e - r Q ) 



Using three solutions r e , r c and r Q to equation f(r) = 0, 



K e = 
K, 



6r e 

(r c - r e )(r c - r Q ) 
6r c 



u (ro-r e )(r c -r ) 

K = A. 

6r Q 



(13) 
(14) 
(15) 



Through tortoise coordinate transformation (TT21 , the radial equation (fTCTj) can be written in the Regge- Wheeler form 

d 2 



dx 2 



AM 2 V(r) 



$ ui (r,) = 4MVf w ,(r,), 



(16) 



which has the form of Schrodinger equation in the quantum mechanics. The incoming or outgoing particles flowing 
between inner horizon r e and outer horizon r c are reflected or transmitted by the potential barrier V(r). It is obvious 
that the evolution of wave solution vf^ of massless scalar field is also determined by potential V(r). 

III. ABSORPTION PROBABILITIES FOR SCALAR EMISSION IN LOW-ENERGY REGIME 



In the above calculation, the radial r is absorbed into wave function ^^(r) to get the Regge- Wheeler equation. 
But in this and next section we need do a reduction run, i.e. using replacement ipuii(r) = rP^iir). So the original 
radial radiation Eq. (|10j) can be rewritten to a new form 



■j- r f(r) — 

dr V dr 



(17) 



It is well known that, in the low-energy regime u — > 0, only the lowest angular momentum will contribute to the cross 
section 



17| . So, in this paper we consider the dominant scalar decay mode S wave (1 = 0), which also is illustrated 
clearly in FigJ5J The analytical solution can be obtained as follows, 

logr log(r-r c ) 



Pui(r) = C 2 + d 



r c r e (rc + r e ) r c (2r 2 - r c r e - r 2 ) 
log(r-r e ) \og(r + r c + r e ) 



7r 2 r e 



7r c r 2 — 2r'f 



r e (—r 2 — r c r e + 2r 2 ) —2r\ 
Using the surface gravities to express above solution, we can get 

log r A log(r — r c ) A log(r — r e ) A log(r + r e + r c ) 



P u i(r) = C 2 + C x 



(18) 



(19) 



r c r e (r c + r e ) 6n c r 2 6n e r 2 6(r e + r c ) 2 K a 

Then we analyse the dynamic behavior of scalar particles propagation near horizons. According to the forenamed 
definition of horizons, we know that one feature of a static spherically symmetric space is that the potential vanishes 
near horizons r e and r c . So we have the relationship /(^)l(r->r e ,r e ) — ^ i n the nearby regions of horizons r e and 
r c outside SdS black hole immediately. Hence, the Regge- Wheeler equation (fTBf can be transformed into a standard 
wave equation form, 

d 2 



dx 2 



- AM 2 uj 2 



**z(r.) =0. 



(20) 



The general boundary condition of greybody factor is that, near the black hole event horizon the solution is purely 
ingoing and near the cosmological horizon the solutions are include incoming and outgoing modes, 



^ NC = B ie -^ r * 



for r ~ r P 



-t>2e , 



fo 



r r ~ r r . 



So near the event horizon r c , the asymptotic solution (1211) yields a low-energy expansion, 



1 AM 



log(r - r e ) log(r c - r) log(r - r ) 



f£c 

Submitting r — >■ r e into solution P w j (1191) , we can get the general solution near event horizon, 



(e) 



Z |r— >r ( 



= C 2 r e + d 



logr e Alog(r-r e ) r e Alog(r c - r e ) r e Alog(r e - r ) 



6r e /t R 



' c' o u ' e^e 

Comparing Eq. (j2~4")) with Eq. (j2"3")l . we can get the expressions of coefficients C\ and C 2 near event horizon, 



C x = -At 



2MA ' 



C 2 = 



(21) 
(22) 

(23) 

(24) 

(25) 
(26) 



So near the cosmological horizon, the low-energy expansion of the asymptotic solution (1221) is obtained by using the 
similar handler routine, 



^(r)c(B 1 +B 2 ) + —(B 2 -B 1 ) 



log(r - r e ) _ log(r c - r) \og(r - r ) 



K e K c Kq 

Similarly, submitting r — > r c into solution P^i (|19p , we can get the general solution near cosmological horizon, 

logr c r c Alog(r c - r e ) Alog(r c -r) r c X log(r c - r a ) 



i> c (r) = rP ul \ r ^r e = C 2 r c + d 



r e r Qr^K e 6r c K c 6r%K 

Hence, comparing Eq. (|2"81) with Eq. (l2"Tl) we obtain the expressions of coefficients G\ and C 2 near cosmological horizon, 



(27) 



(28) 



3 1 

C x = -iwr c (B 2 - Bx)——, 
2 K 1 MA. 

C 2 = -(Bx+B 2 ) + 0{uj). 



(29) 
(30) 



Using the solutions near event horizon and cosmological horizon, we then stretch and match them in the intermediate 
region and obtain the key absorption coefficients. Distinctly, the final relationship of Ax, Bx, B 2 can be obtained 
according to Eqs.®,®,®,®, 

Ax 



Bx+B 2 
Ax 



B 2 — Bx 

The corresponding reflection and absorption probabilities are given naturally as follows, 

,B2 l2 (r 2 c -rlf 



\W 



r 2\2 ' 



„2„2 

1-I^| 2 



At v 
^' e ' c 



{r 2 c +rl) 2 ' 



(31) 
(32) 

(33) 
(34) 
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FIG. 4: absorption probabilities \A\ 2 versus cosmological constant A for unit mass M = 1 in low-energy regime. 



36, 



37 



52j . Here, we present the solution in alternative 



This result is in accordance with that of two-horizon spacetimes 
method. In order to get the variation of absorption probability with cosmological constant, we use expression Eq.® 
to replace r c and r e in above Eq.(|34j). Finally, the analytic expression of absorption probability is obtained as follows 



2 _ 12 + 8 cos Atjq - 8 cos 2( - 8^3 sin 2( 

_ 18 + 4cos4r? - 8cos2C + cos4C- 873sin2C- y3sin4C' 

where the new variables 770, C arc defined by r/Q = l/3arccos (3MVA), C = 1/3 arcsin(3Mv / A) respectively. One 
interesting fact of low energy regime is that this kind of scattering problem varies only with the whole value of M-\/A 
which has a dimension of GeV 2 coincided with the reciprocal quantity of Newton gravity constant. The form of 
the absorption probabilities for low-energy regime is shown in Fig. It illustrates that absorption probabilities 
\A\ 2 increases monotonically with cosmological constant. The maximum cosmological constant 0.11 is corresponding 
maximum probability |*4|ff Q:! . = 0.986701. This point also matches the potential V(r) very well (please see Fig. 
[3]). Because that potential becomes lower with increasing cosmological constant, the action of lowering potential 
on the particles is weak and the absorption probabilities are naturally larger. On the other hand, it is shown that 
the absorption probability is vanished when cosmological constant disappears. This point is also justified by papers 



35 37 



52j. Referring to the cross section a s ([T]) and transmission coefficients (T s )i ((2|) in 4D Schwarzschild space, 
we can see that when the particles mass is zero and the energy to is more less than the mass of black hole 2M, the 
cross section reduces to <j s ~ IQirM 2 . This result shows that, in the low-energy regime the absorption cross section 
matches accurately with the area of Schwarzschild black hole. This point is expanded to arbitrary dimensions by Das 
et.al 17j latterly. From the Eqs. ([1]) and ([2]), we can also see that the transmission coefficients (T s )i vanish in the 
low-energy limit w — > for all values of angular momentum / including s wave, which justify the induction result 
of zero cosmological constant in SdS black hole. Unlike the singular horizon in Schwarzschild space, there are two 
different horizons (r e , r c ) in Schwarzschild-de Sitter space and the second horizon creates a finite universe. So for 
the low-energy particles, which has a infinite wavelength, can be localized and has a finite probability to traverse the 



finite distance between the two horizons 13 
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FIG. 5: effective potentials 1/B 2 (r) of high-energy particles with various cosmological constants and unit mass M = 1. 



IV. GREYBODY FACTOR OF SCALAR EMISSION IN HIGH-ENERGY REGIME 



When the particles pro pag ate in the high-energy regime, the scalar greybody factor can be assumed its geometric 



optics limit value 
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16, 
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68H7Gl| . The momentum of high-energy particle is p^ 



dx* 

dc ■ 



where C is 



an arbitrary scalar parameter for the affine parameter. So the first and second motion equations about t and <p are 
obtained naturally in the following equations, 

2 dp 



a . dt 



And the third motion equation can be given by the normalization relation of particles g^p^p" = 0, 



L 2 



UJ 2-11 



2M A ; 

r 3 



(36) 
(37) 

(38) 



Then the high-energy particle is subjected Eqs. ((36|) . (1371) and (|38|) near SdS black hole. Combining Eq. (|36|) with Eq. 
\, the trajectory equation can be gotten as, 

2 



1 dr 



K ~, I = ^ 1- 



1 



2M A 



(39) 



v r 2 dtp J \LJ r 1 \ r 3 

Here, we define two new parameters, one is the impact parameter b — L/u>, which means the effective sighting range, 
the other is the high-energy particle's effective potential 1/B 2 (r), which is shown in FigJS] The parameter B{r) is 

'1/2 



2M A 9 

B(r) = r 1 r 2 

r 3 



According to Eq. (|40|) . particles orbital equation (139)) can be rewritten as 

2 



1 dr 
r 2 df 



1 



1 



b 2 B 2 (r)' 



(40) 



(41) 
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FIG. 6: greybody factor versus cosmological constant A for unit mass M = 1 in high-energy regime. 



which is also called Binet equation in the subject of celestial bodies motion orbits. The classically accessible regime 
is bMax = B(r)Min, which means the closest orbit of high-energy particle before absorption by black hole. The 
absorption area of the black hole with high energies is 

-l 



a h = Tib 



Ma 



1 



27 M 2 » ' ■ (42) 
which is shown in Fig|6l Like the regime of low energy, the greybody factor ah is also increasing monotonically 
with increase in cosmological constant A in high-energy regime. When the cosmological constant vanishes, the above 
expression reduces to the Schwarzschild case ah = 27nM 2 = 27nr 2 /4: 



18, 



19 



30 
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CONCLUSION 



In this paper, we have studied the scalar particles radiations outside of Schwarzschild-de Sitter black hole in low- 
energy regime and high-energy regime. Finally, we have successfully obtained the analytic absorption probability 
\A\ 2 in low-energy regime and the greybody factor ah in high-energy regime, respectively. We find the cosmological 
constant influences intensely the cross section and greybody not only in the low-energy zone but also in the high-energy 
zone. 

Firstly, we have a look at the case of low-energy limit. It is well known that the cross section of Schwarzschild 
black hole vanishes in the low-energy and low angular momentum limit 13j. However, in the SdS case the whole 



space is deformed due to the existence of cosmological constant that changes naturally the scatting of scalar field. 
For example, the improved potential V(r) (jlip . which is illustrated in Figs.©, decreases with increasing A. So 
the absorbent capacity, which expresses mathematically by absorption probabilities \A\ 2 , enhances for a low potential 
V(r)i ow . This point also explain essentially the tendency profile of Fig.Q, which shows absorption probabilities \A\ 2 
increases with increasing in cosmological constant. 

Secondly, we make some discussions on the high-energy limit. Due to the special physical characteristics of high- 
energy scalar particles, it is appropriate to assume the geometric optics limit value. Likewise, the greybody factor ah 



11 



also is under the influence of A. Similarly, the stronger absorption capacity for increase in cosmological constant A 
determines the same monotonicity for the greybody factor 07,. in low-energy regime in Fig. ([6]). The only difference 
between low-energy and high-energy regimes is that, for the large cosmological constant the ah increases more rapidly. 
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